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Abstract. We review the general interplay between Nuclear Physics and neutrino-nucleus cross sections at intermediate and 
high energies. The effects of different reaction mechanisms over the neutrino observables are illustrated with examples in 
calculations using several nuclear models and ingredients. 
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As motivation for the more specific workshop sessions, in this talk we introduce the general formalism of neutrino 
scattering from nuclei and define the observables of interest for nuclear physics at the energy regime of interest. Using 
different nuclear models we present with examples the theoretical ingredients of relevance for neutrino reactions: 
long Range nuclear correlations (RPA), final state interactions (FSI), finite-size effects, Coulomb corrections, and 
relativistic effects. Theoretical results will be shown for charge-changing quasielastic neutrino scattering. We skip here 
the discussion on neutral current scattering, Delta excitation and coherent pion production, which will be summarized 
in the M.B. Barbara and E. Hernandez papers in these proceedings. We present results for kinematics going from 
low to high energy and for different kind of observables: response functions, inclusive cross sections, integrated cross 
sections, angular distributions, polarization observables, etc. Nuclear models for which we will show results are Local 
Fermi Gas (LFG), Relativistic Fermi Gas (RFG), Semi-relativistic Shell Model (SR-SM), and Super-Scaling Analysis 
(SuSA) model. We also skip the discussion on the Relativistic Mean Field, which will be summarized in the talk by 
J.M. Udias in these same proceedings. Some particular topics that we briefly discuss are theoretical uncertainties on the 
ratios of interest for experiments on atmospheric neutrinos, nuclear effects on lepton polarization, and reconstruction 
of neutrino cross section from electron scattering data. 



Formalism 

We consider as example the CC neutrino reaction V/ +A — * l~ +B, for initial and final K 1 ^ lepton momenta, with 
energies e and e', respectively. We introduce also the momentum transfer Qj^ = — K ,f * = P ,fi —P^.lt is convenient 
to define the following adimensional variables A = K = ^j-, and x = K 2 — X 2 . 

The inclusive cross section for the inclusive (v/,Z~) reaction where only the final lepton is detected, can be written 

as 

da G 2 cos 2 d c , , , ? 9 _ 7 
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Where the coupling constant G = 1 .1664 x 10 5 GeV , the Cabibbo angle C = 0.974, and the generalized scattering 
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angle is defined by tan 2 f = -. J®J_ 2 ■ The interesting nuclear information is contained in the structure function 



&l = VccRcc + 2V cl Rcl + V LL R LL + V T R T + 2V T ,R r . 

The kinematic factors Vk, coming from the leptonic tensor, are defined by |Q]] 
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V cc = 1-S 2 tan 2 -, V LL = ^+ 1 + — +p5 2 U 2 tan 2 -, y CL = - + — tan 2 - 
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where the following adimensional variables are introduced: 8 = , p = and p' = — 
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Finally the following five nuclear weak response functions appear 

R CC = W 00 , R U = W 33 , RcL = -\{W m + W 30 ) 

R T =W n +W 22 R T , = - l -(W n -W 2X ) 

as linear combinations of the weak CC hadronic tensor. 

W» v ^(Ef-Ei-coXflJ^qWifirWli) 

fi 

where the matrix elements of the hadronic current are taken between the initial and final hadronic states in the nuclear 
transition |/) — * |/). The hadronic current for a single nucleon is of the form V — A 



J PL = u p {p') 



F l (Q 2 )^+F2(Q 2 )i^ v ^--G A (Q 2 )^Y5-G P (Q 2 )^-Y5 



The Relativistic Fermi Gas (RFG) 
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As example of the above general formalism in a many-nucleon system, we present here the nuclear response 
functions in the RFG, where the initial and final nucleons are described as (relativistic) plane waves. For the inclusive 
(V/,/~) reaction we can write the five response functions as 

RK=NA UKfRFG(Y) = G K fRFG(Y), K = CC ,CL,LL,T ,T' , (2) 

Here N is the neutron number, and we have defined An = — with 77 f = kr Imv is the Fermi momentum in units 



of the nucleon mass, and %p = y 1 + t]p — 1 is the Fermi kinetic energy measured in the same units All the response 
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functions are proportional to the RFG Scaling function /rfc^Y) = 1(1 — V^ 2 ) ^ ( 1 — W 2 )> tnat depends only on the 



scaling variable 

1 A-T 
W= -7*=-, ; (3) 

V?f v /(1+A)t+k-Vt(1 + t) 
Explicit expressions can be obtained for the single-nucleon responses Uk id. 



The Local Fermi Gas (LFG) 



As an extension of the constant-density Fermi gas, the LFG is an easy model to include non-trivial nuclear effects 
J2l . It is based on the local density Approximation (LDA), with a local Fermi momentum is kp (r) = (3 lt 2 p (r) ) 1 / 3 . The 




FIGURE 1. Left: LFG predictions for neutrino reactions from 12 C. Left panels: (Vn Right panels: (e, V e ). Pauli+Q are the 
LFG results with correct energy balance. RPA are the results with RPA and Coulomb corrections. Top panels: total cross sections 
multiplied by the neutrino fluxes of the bottom. Middle panels: differential cross section for fixed neutrino energy. Right top: v f 
and V e differential cross section for fixed momentum transfer; right bottom: effect of the FSI over the neutrino cross section, with 
and without RPA. 



TABLE 1. Flux averaged neutrino cross section 
from 12 C in 1CT 40 cm . The experimental data are 
fromLSND HH. 





Pauli+Q 


RPA 


LSND 


) 

(e, v e ) 


20.7 
0.19 


11.9 
0.14 


10.6±0.3±1.8 
0.15±0.01±0.01 



responses are averaged over the nuclear interior, weighted by the nucleon density p(r). In this model one can easily 
include relevant nuclear effects such as the correct energy balance RPA nuclear correlations, Coulomb distortion, and 
FSI effects. 

While the LFG correctly takes into account Pauli-blocking effects, it gives a wrong energy balance for the nuclear 
excitations. The energy balance can be corrected by the minimum nuclear excitation energy gap Q = M(Xf) — M(Xj), 
instead of the usual LFG value Q.LFG{ r ) = Ep (r) — Ep(r). That is equivalent to replacing 0) — ► (0 — [Q — Qlf g ('")]■ 

The RPA series of Fig. 1 is solved by using a ph-ph interaction of Landau-Migdal, with parameters fitted to 
electromagnetic nuclear properties and transitions [3]. The RPA series can be generalized by including A excitations 
in the medium, and ph-Ah, A-hAh effective interactions. The sum of the RPA series is equivalent to a renormalization 
of the axial and vector parts of the weak hadronic tensor in the medium 

Coulomb corrections can be included in the LFG by introducing the Coulomb self-energy of the final lepton 
T.C = 2e'Vc{r) (where Vc(r) is the nuclear Coulomb potential) in the charged lepton propagator, k2 ^_ ^ - ^ +iE - 

A new local energy-momentum relation for the final lepton is then obtained, that is used in the LFG calculation. This 
procedure is equivalent to the modified effective momentum approximation. 

LFG results for 12 C(Vjj,/i~) and 12 C(v e ,e~) are presented in Fig 1, and are compared with experimental data 
10, Eh in Table 1. The LFG with correct energy balance (Pauli+Q), over-estimates the data, while a good agreement is 
achieved once RPA and Coulomb corrections are included. 

Final State Interaction (FSI) can be taken into account in the LFG model by using a renormalized nucleon propagator 
in the medium Gfsi(p) — o_wi_w j , where £(/?) is the nucleon self-energy in the medium. A good approximations 
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FIGURE 2. Ratio of inclusive cross sections for argon. The 68% confidence level band is displayed for the full LFG model 
including RPA, Coulomb and FSI. We also show results with the bare Fermi gas model without nuclear corrections. 

for intermediate energies is to take Im£/, ~ for hole states. The nucleon self energy can be computed by diagrammatic 
techniques [6]. As shown in Fig. 1, the FSI importantly changes the shape of the differential cross section. The main 
effect is is an enhancement of the high energy tail and a reduction of the cross section at the peak region. In the same 
figure we also show that the RPA corrections are less important in presence of FSI. 

Theoretical uncertainties in the LFG model have been computed assuming central values and errors of the model 
input parameters, and including 10% uncertainties in both the real part of the nucleon self-energy and densities |7[]. 
A Montecarlo simulation is then performed by generating sets of input parameters using Gaussian distributions. After 
computing the different observables, one obtains the distribution of the observable values, and the theoretical errors are 
identified by discarding the highest and lowest 16% of the obtained values, keeping a 68% confidence level interval. 
Uncertainties on the integrated cross sections are of the order of 10-15%, which turn out to be similar to those assumed 

for the input parameters. As shown in Fig. 3, theoretical errors cancel partially out in the ratio — of interest 

for experiments on atmospheric neutrinos. 



Polarization observables 

The study of final % polarization in (v T , i) reactions |8] is of interest for Vu — ► V T oscillation experiments. The z 
decay particle distribution depend on the z spin direction, and thus theoretical information on the z polarization will 
be valuable OJ. Information about Z polarization is also needed in Vu — ► V e oscillation experiments to disentangle 
(v e ,e) events from background electron productions following the — > V T oscillation [ 10]. 

The polarized differential cross section in (V/,Z) reactions when the final lepton polarization is measured in the 
direction s can be written as 



E(3) = 



d 2 o 
dCl'dE' 



1 



where Eo is the unpolarized cross section. The lepton polarization vector components P/ (longitudinal component, in 
the direction of the final lepton), and P, (transverse component in the scattering plane) can be obtained as asymmetries 



E(S) +£(-*)' 



Results for the polarization observables with the LFG model are shown in Fig. 4, with the differential cross section, 
the two T polarization components, and the total polarization and angle. The nuclear RPA and FSI effects are of small 
importance at the quasielastic peak region because of partial cancellation when computing the asymmetries. Their 
importance increases in the tail of the cross section. 
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FIGURE 3. Polarization observables in the LFG (dotted), adding RPA correlations (dashed), and adding FSI (solid). 



Super- Scaling Analysis (SuSA) 

In the RFG, Eqs. (12131 1. all the response functions can be factorized in a single-nucleon function times the superscal- 
ing function /rfg* which only depends on the scaling variable \j/ and is the same for all nuclei. We then say that the 
RFG superscales, namely /rfg does not explicitly depends upon the momentum transfer q (scaling of the first kind) 
nor the Fermi momentum kf (scaling of the second kind). Moreover the scaling function is the same for electron and 
neutrino scattering, hence the neutrino and electron cross sections are related just by a single-nucleon factor. 

These ideas have been extended to extract an "experimental" scaling function from (e,e') data by computing the 
ratio 

da 



fexp{V) = 



dSl'ds' 



exp 



GyLott { v lGl + vtGj) 



where is the scaling variable shifted to account for an energy binding parameter E s . An extensive analysis of (e,e') 
data has been performed in the quasielastic peak fiUl . The parameters kp y E s are fitted to the data to minimize 
the differences between scaling functions for different kinematics and nuclei. When this analysis is performed on 
the experimentally available longitudinal, fi = and transverse, fj — response functions, one finds that the 
L response function superscales quite well, while scaling is broken in the T response above the peak, due to non- 
quasielastic processes (pion production, A excitation, MEC, etc). The experimental scaling function is shown in Fig. [6] 
below 

Having succeeded in representing the (e,e r ) quasielastic data by means of an universal superscaling function, 
the SuSA procedure can be reversed and give predictions for neutrino reactions lfl2tl . In fact, starting with the 
experimental (e.e 1 ) scaling function one can use the RFG equations (0 to compute the (V/,/~) response functions 
with the substitution /rfg ( W) — * fexp ( W) ■ The assumption underlying this procedure is that the superscaling function 
associated to the 5 weak responses is the same and is equal to the longitudinal electromagnetic response function. 
Calculations aiming to to justify theoretically the validity of SuSA are presented in the next section. 
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FIGURE 4. Scaling properties of the SR shell model. Left: scaling of the first kind. Curves for q = 0.5, 0.7, 1,1.3,1.5 GeV 
collapse into one. Right: scaling of the second kind. Curves for 12 C, 16 and 40 Ca collapse into one 

The semirelativistic shell model (SRSM) 

In order to study the scaling properties in realistic models one needs to go beyond the simplicity of the Fermi gas, 
where scaling holds by construction. In the continuum shell model (CSM), i.e., nucleons in a mean field, distortion 
of the ejected nucleon is present and a general proof of scaling cannot be provided. Thus, at least within the context 
of the CSM, we shall be able to check the consistency of the SuSA approach and quantify the degree to which scale 
breaking effects are expected to enter. 

In the CSM the initial state |/) is described by a Slater determinant with all shells occupied, while in the impulse 
approximation the final states are particle-hole excitations coupled to total angular momentum \f) — \{ph )J). The 
single hole wave function is then written as \h) = \£hhjh) while the single particle wave function \p) — \e p l p j p ). The 
radial functions are obtained by solving the Schrodinger equation with a Woods-Saxon potential. 

At the ongoing and next generation neutrino experiments the neutrino beam energies increase to the GeV level and, 
typically, large energies and momenta are transferred to the nucleus. For these kinematics, relativity is important. The 
relevant relativistic corrections can be easily implemented in the CSM through the semi-relativistic approach (SR) 
[Q]]. It is based on a new expansion of the relativistic single-nucleon current j^(p',p) = u(p')r fl (Q)u(p) in powers 
of the initial nucleon momentum, r\ = p/m^, to first order 0(r\). We do not expand in p' /ni/j, hence q and CO can be 
arbitrarily large. Second one must use relativistic kinematics. The energy transfer in the CSM is the difference between 
the (non-relativistic) single-particle energies of particle and hole co = E p — £/,. The relativistic kinematics are taken into 
account by the substitution e p — > £ p (l + e p /2niN) as the eigenvalue of the Schrodinger equation for the particle in the 
continuum. 

A test of the SR approach was performed in fljj] by comparing the exact RFG results with those of the SR Fermi gas 
model, where the above SR current and relativistic kinematics have been implemented. 

An study of superscaling properties of the SRSM, for q > 0.5 GeV/c, is performed in Fig. [4] Scaling of both first 
and second kind are achieved in the shell model: all the curves collapse into one with small deviations for q = 0.5 
GeV. Since both kind of scaling are found, we conclude that superscaling occurs within our model. 

The SRSM with Woods-Saxon potential fails to reproduce the experimental scaling function, giving results similar 




FIGURE 5. Scaling functions computed in the SR shell model compared with the relativistic mean field model (RMF). Dotted: 
Woods-Saxon potential. Solid: DEB+D potential. Dashed: RMF. 



to the RFG, with the exception of small tails for low and high energy transfer. The next step lfl3ll is to improve the 
relativistic description of the ejected nucleon. The FSI is modified by using a Dirac-equation based potential plus 
Darwin term (DEB+D) in the final state, instead of the usual Woods-Saxon potential. The trick is to rewrite the 
Dirac equation as a second-order equation for the upper component y/ up (f). The Darwin term is then defined by: 
Vupir) — K(r,E)<j)(r), where the function <j>(r) verifies the Schrodinger equation 



' V 2 + U DEB (r,E) 



2niN 



Both the DEB potential Udeb^E) and Darwin term K(r,E) are energy-dependent functions. The scaling functions 
computed using the DEB+D potential are compared to the exact RMF results in Fig. [5] Both models give similar 
results, and clearly different to the SR model based on a Woods-Saxon potential. 

Scaling of the first kind of the SRSM with the DEB+D potential is shown in Fig. |6j where we also compare with 
the experimental data for fa. The scaling is surprisingly good, taking into account that the momentum transfer ranges 
from 0.5 to 1.5 GeV. A q-dependent energy shift E s (q) is needed to bring all the curves to collapse. The q dependence 
of E s {q) is linear. 

A test of the SuSA reconstruction of the (v^ju - ) cross section from the (e,e') one is shown in Fig. [6] There the 
neutrino cross section obtained by direct computation is compared with the SuSA reconstruction from the computed 
(e,e') longitudinal scaling function. Both results are quite similar for scattering angle above 30°, corresponding to 
intermediate to large momentum transfer. For very low scattering angle there are small differences between both 
procedures due to the small momentum transfers involved. 



Concluding remarks 

Neutrino interactions importance for nuclear physics has been illustrated with examples from a selection of neutrino 
reactions on nuclei using different approaches. The connection of electron and neutrino cross sections has also been 
analyzed within the superscaling approach. Since neutrino cross sections incorporate a richer information on nuclear 




FIGURE 6. Left: Scaling of 1st kind with DEB+D potential compared with experimental data, q = 0.5, 0.7 1.0, 1.3 and 1.5 
GEV/c. Right: Test of SuSA in the SR shell model for the 12 C(Vjj,/Z _ ) reaction with neutrino energy e = 1 GeV. Dotted: Woods- 
Saxon potential. Solid: DEB+D potential. Dashed: SuSA reconstruction from the computed (e,e') scaling function. 



structure and interactions than electron cross sections, the availability of neutrino-nucleus observables of different 
kinds will be valuable for the development of more precise nuclear models and nuclear interaction theories. 
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